Larmor precession and Debye relaxation of single-domain magnetic nanoparticles 
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The numerous phenomenological equations used in the study of the behaviour of single-domain 
magnetic nanoparticles are described and some issues clarified by means of qualitative comparison. 
To enable a quantitative application of the model based on the Debye (exponential) relaxation and 
the torque driving the Larmor precession, we present analytical solutions for the steady states in 
presence of linearly and circularly polarized ac magnetic fields. The power loss is calculated for both 
cases and the results compared with the predictions of the Landau-Lifshitz-Gilbert equation. Using 
the exact analytical solutions, we can confirm the insight that underlies Rosensweig's introduction 
of the "chord" susceptibility for an approximate calculation of the losses. We also find that this 
approximation provides satisfactory numerical accuracy up to magnetic fields that would generate 
3/4 of the saturation magnetization if applied constantly. 

PACS numbers: 47.65.Cb, 75.30.Cr, 75.75.Jn 
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I. INTRODUCTION 

Enduring interest in colloidal dispersions of ferro- and 
ferrimagnetic nanoparticles is fuelled by their medical ap- 
plications. Colloids of small particles of iron oxide are 
used in magnetic resonance imaging (MRI)i as contrast 
agents and in hyperthermia treatment as heat-generating 
media. 2 In both cases, the nanoparticles are subjected to 
a time-dependent magnetic field, but the frequencies and 
magnetic field intensities are quite different. In the case 
of MRI they are determined by the resonance frequency 
of proton, which is 63 MHz in a field commonly used, 1.5 
Tesla. In hyperthermia, medical considerations require 
frequencies of the order 10 5 Hz and magnetic fields of 
10~ 2 Tesla. The frequency of ferrimagnetic resonance in 
magnetite is higher than 1 GHz at zero external field, 3 so 
that hyperthermia takes place at frequencies way below 
resonance, irrespective of the applied field. 

The theoretical background of the processes taking 
place in these applications is colorful. Various approaches 
are being pursued with little attention for attempts to 
compare and evaluate the available methods. Berger et 
alA have undertaken to compare eight phenomenologi- 
cal equations of motion for magnetic moments of ferro- 
or ferrimagnetic nanoparticles dispersed in a nonmag- 
netic medium. They found that the calculated line- 
shapes in the narrow-linewidth limit, pertaining to high 
temperatures, are of the Lorentzian form, irrespective 
of the equation of motion used, but in the case of 
broad resonance lines significant variations are observed. 
Low-temperature superparamagnetic resonance spectra 
proved to enable a distinction between different equa- 
tions of motion, of which only the Landau-Lifshitz equa- 
tion provided a satisfactory fit. 

The purpose of the work reported here is a theoreti- 
cal comparison of two equations of motion, the Landau 
-Lifshitz-Gilbert equation and the modified Bloch equa- 
tion (in the nomenclature of Berger et al£). Our focus 
is on hyperthermia, rather than magnetic resonance. We 
are seeking analytical solutions to the equations of mo- 



tion in the absence of a strong static magnetic field, which 
give, in addition to the specific losses, some insight into 
the underlying physical processes. 



II. EQUATIONS OF MOTION AND 
SIMPLIFICATIONS 

The behavior of single-domain ferro- or ferrimagnetic 
nanoparticles in an external magnetic field has much in 
common with that of atomic or nuclear magnetic mo- 
ments. The torque on a magnetic moment fi, 



T = /xx B 



(1) 



determines the equation of motion of the angular mo- 
ment, dL/dt = T. The gyromagnetic relation, fi = 
enables a closed equation for fj, which, applied to the 
magnetic moment of unit volume, provides the equation 
of motion of the magnetization, 



dM/dt = 7M x B. 



(2) 



Here 7 is the gyromagnetic ratio. The magnetization 
of materials we have in mind in this work is due to the 
spin of electrons, accordingly 7 = —1.76 x 10 11 Am 2 /Js. 

The vector product in eq.([2| implies that any change 
of the magnetization is perpendicular to M , that is, the 
modulus of M remains constant. Also, if B is constant, 
d(M ■ B)/dt — 0, that is, the angle between the two 
vectors is constant. The only motion satisfying these 
conditions is a precession of M around B. The angular 
velocity of the magnetization in this Larmor precession 
is u)l = 7 M x B I My , where Mj_ is the projection of M 
on the plane perpendicular to B. The Larmor frequency 
is defined as a positive quantity, wi = |7|B. To reduce 
the potential energy, U = —M ■ B, the contribution to 
dM/dt due to relaxation must have a component parallel 
with B. Landau and Lifshita^ have chosen a "damping 



2 



term", which evidently achieves this, being proportional 
to M x (M x B) = (M ■ B)M - M 2 B. The Landau- 
Lifshitz equation of motion is then 



dM 
~~dT 



= j[MxB + aM~ 1 Mx(MxB)]. (3) 



The coefficient a goes under the name of "the dimen- 
sionless damping coefficient", which is something of a 
misnomer, because addition of the "damping term" evi- 
dently enhances the motion of M: 



dM 



dt 



= j\M x B\(l + a z ) 



2x1/2 



(4) 



Gilbert's approach 6 is closer to the notion of friction, 
as it subtracts from the Larmor torque a torque propor- 
tional to \dM/dt | . This is reminiscent of friction in linear 
motion, where a force opposite to \dr/dt\ is introduced 
into the equation of motion. The analogy allows a deriva- 
tion of the Gilbert equation, 



dM 
~dT 



= 7[M x B - rifj, M x dM/dt], 



(5) 



by adding a Rayleigh dissipation function to the La- 
grangian which describes the Larmor precession of the 
magnetic moment. 

A comparison of the Landau-Lifshitz and Gilbert equa- 
tions reveals that dM/dt (i) is perpendicular to M in 
both equations and (ii) in the former it consists of two 
mutually perpendicular terms while in the latter this is 
not the case. Observation (ii) implies that decomposition 
of the damping term in the Gilbert equation into compo- 
nents parallel and perpendicular to that of the Landau- 
Lifshitz equation will offer a direct comparison of the 
two. In fact, the Gilbert equation itself provides a de- 
composition into components which delivers the desired 
transformation. Multiplying both sides of eq.© by M 
and taking (i) in account yield 



M x = j[M x (M xB) + n^M 2 dM/dt]. (6) 

Substituting this result in the last term of the Gilbert 
equation and rearranging terms lead to the Landau- 
Lifshitz- Gilbert (LLG) equation, 



dM 
~~dT 



7(1 + a 2 )~ 1 [M x B — aM~ 1 M x (M x B)}, (7) 



where a = "fn/ioM. Clearly, for a -C 1 the Landau- 
Lifshitz equation is a good approximation, but for the 
general case the (1 + a 2 )^ 1 factor is essential to elimi- 
nate the non-physical implications of the Landau-Lifshitz 
equation pointed out by Kikuchi 7 and Gilbert. 6 Also, due 
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FIG. 1: (Color online) Time dependence of M(t)/M s is de- 
rived from the solution of LLG equation © (curve (a)) with 
condition oiujl = 1- For comparison an (1 — exp -1 ) type ex- 
ponential relaxation is also given (curve (b)). Both saturate 
to the same Ms but the path is different. 



to this factor, Gilbert's damping term reduces the motion 
of M: 



dM 



dt 



7|M x B\ (I + a 2 



-1/2 



(8) 



Strictly speaking, the effect of the Landau-Lifshitz and 
Gilbert damping coefficients cannot be described with a 
relaxation time, because the relaxation they stand for 
is not exponential. If B is a constant field, pointing 
in the z direction, the solution of the LLG equation is 
M z = Mtanh(ao)Lt), with Co L = uj l /(1 + a 2 )~ 1/2 . If 
the magnetization is aligned perpendicular to B at the 
outset, switching into the energetically favourable state 
proceeds as given in Figfl] curve (a). For comparison, 
curve (b) shows the path of exponential relaxation with 
the same initial state. The similarity may suggest the 
definition of 



tllg 



auj L 



(I + a 2 ) 1 / 2 



(9) 



as the relaxation time of the LLG dynamics. However, 
its dependence on B would make tllg an awkward pa- 
rameter when the magnetic field is time-dependent. Fur- 
thermore, a constant value for the product tllg&l, is 
counter-intuitive as different values of the ratio of the 
relaxation time and the Larmor precession time are ex- 
pected to lead to qualitatively different behavior. 

Shliomis 8 has suggested that under well-defined condi- 
tions the equation of exponential relaxation, 



dM 
~~dT 



M M P . 



(10) 
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should suffice to describe the behavior of a colloid of mag- 
netic nanoparticles. Here, M is the average magnetiza- 
tion of the particles, 



M e . 



M S C 



( (JLoHMgV 

\ kT 



en, 



(11) 



V is the particle volume, en = H/H is the unite vector 
pointing along H and M$ is the saturation magnetiza- 
tion of the colloid Ms — <pMd, where <fi is the volume frac- 
tion and Md is the magnetization of the single-domain 
magnetic nanoparticle. 

The Langevin function, C(x) — coth(a;) — gives 
the magnitude of the magnetization in thermal equilib- 
rium. Note that M eq must be an ensemble average and 
consequently so is M. In this respect, in the context of 
superparamagnetic resonance or hyperthermia , eq. (|10[) 
is more expedient than the LLG equation. In the latter 
case, having found the possible solutions of the equation 
of motion, one has to face the issue of the appropriate 
weighted average of the associated energy losses. 

Clearly, the parameter r in eq. ljlOp is a proper relax- 
ation time. In a stationary field, where M eq is also con- 
stant, the solution of this equation for the component of 
M along H is 



(M(t) - M eq ) = (M(0) - M eq ) exp(-t/r). 



(12) 



In reference^, Shliomis names eq. (|10[) the Debye re- 
laxation equation. Yet it is often reduced from what is 
generally called the Shliomis relaxation equation. 



dM 
~~dT 



M (V ■ v) - oj x M : 



M - M„ 



(13) 



Devoid of the possibility of Larmor precession, the De- 
bye relaxation equation, eq. (|10|) . can be looked upon 
as a truncated and simplified version of the Bloch- 
Bloembergen equation. The latter does include the gyro- 
magnetic torque, [IqjM x H. The equation of motion, 
which Bloch has set up for the interpretation of nuclear 
magnetic resonance experiments^ and Bloembergen ap- 
plied to ferromagnetic resonance^ was specifically aimed 
at the configuration of such experiments. The z axis is 
aligned along the static magnetic field, which provides 
the axis of the Larmor precession, and enables the defi- 
nition of two distinct relaxation times: t\ for the z com- 
ponent of the magnetization vector and T2 for its compo- 
nents in the xy plane. As the Debye relaxation equation 
is not meant to be limited to ferrofluids in a dominant 
static magnetic field, the simplification to a single relax- 
ation time seems inevitable. 

Usually, the equilibrium magnetization in the Bloch- 
Bloembergen equation is assumed to be generated by the 
static fields. The effect of the much weaker microwave 
field on the equilibrium magnetization is neglected. In 
what Bcrger et alX call the modified Bloch equation, the 
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FIG. 2: Color online. The variation of exact equilibrium and 
'chord'-susceptibility based magnetizations at field amplitude 
of H = 50 kA/m. 



equilibrium magnetization is defined in terms of the to- 
tal effective magnetic field, which includes the microwave 
field as well as the internal fields (demagnetizing field 
and crystalline anisotropic field). The latter are rea- 
sonable, internal fields are often included in the Bloch- 
Bloembergen equation as well. However, the legitimacy 
of a time-dependent field in the same place implies the 
assumption that the response of the magnetization to a 
rapidly changing field can be described in the same way 
as in the case of a stationary field. 

Recently, Cantillon-Murphy et alr^ have published a 
thorough analysis of the implications of the Shliomis 
relaxation equation for the relaxation process and en- 
ergy dissipation, in superparamagnetic colloids contain- 
ing magnetite nanoparticles, under the conditions of 
magnetic resonance imaging (MRI) and hyperthermia 
treatment of cancer patients^ They justify the neglect 
of v and u>, in eq. lfTB")) so that in fact the calculations are 
based on the Debye relaxation equation, eq. tfTU)) . Apart 
from the approximation inherent in the exclusion of 
the gyromagnetic torque they also applied Rosensweig's 
chord susceptibility^ 



Mg r ( mH Q M d V 
Xch ~ H I kT 



(14) 



instead of the Langevin function appearing in eq. (|ll[) . 

The use of the chord susceptibility in the present con- 
text is illustrated in FigO As the alternating applied 
field is increasing towards its maximum at Hq, eq. ljlip 
requires the equilibrium magnetization to follow the con- 
cave curve denoted by C, whereas the approximate equi- 
librium magnetization given by the chord susceptibil- 
ity increases along the straight line between M eq (0) at 
M eq (Ho). The equation of that line is 
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FIG. 3: (Color online) The external magnetic field depen- 
dence of the behavior of equilibrium and approximate equi- 
librium magnetizations for a single H cycle. The equilib- 
rium magnetizations defined with Rosensweig's susceptibility 
and Langevin function are denoted by R (dotdashed) and 
C (dashed) respectively. At weak fields the two curves are 
identical Fig|3][a). Increasing field strenght only the extrema 
remain identical. 
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FIG. 4: (Color online) Equilibrium magnetization (£) and ap- 
proximate equilibrium magnetization calculated for t — 7T/8 
and t — T, plotted against the amplitude Ho of the oscillating 
magnetic field. At t = T the two are identical by definition 
of Xch, but at t = 7T/8 the relative difference between the 
(£) curve (down triangles) and the (R) curve (up triangles) 
increases with Ho, reaching 20% at 100 kA/m. 



M s -£ ^ fcT ) e ff , (15) 



hence the definition of the chord susceptibility given by 
eq-dHJ). 

Equation (jTUJ) implies that in a field H(t) = 
H cos(27rf/T) the relaxation pulls the magnetization to- 
wards the equilibrium magnetization, which is oscillat- 
ing in time. Figure [3] shows this equilibrium magnetiza- 
tion determined by the Langevin function (£) and the 
Rosensweig chord susceptibility (R). By definition of the 
latter, the two curves meet at H = and H = Hq. In 
the limit of strong external field (fioHoMaV/kT 3> 1), 
Figj3jd), the C curve is flipping between the extrema, 
while the R curve has sinusoidal characteristic. The 
magnetization does not exceed the saturation magneti- 
zation Ms, even for field Hq > Ms- In the weak field 
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limit, Fig[3](a), where hyperthermia is applied, there is 
no difference. Figures (3Jb) and[3Jc) show the transition 
between the two behaviors. It is clear that when the 
field reaches a value where M eq /Ms > 0.8 (Figs. [3Jc) 
and(3Jd)), the magnetization calculated with the chord 
susceptibility is substantially deviates from the correct 
function. 

The dependence of the discrepancy between the two 
curves on the amplitude of the ac field is shown in FigJU 
In the low-field limit there is no discrepancy, because 
Xch = X an d M eq = \H , but it is visible at 20 kA/m 
and increases with increasing field. Beyond Ho = 100 
kA/m the relative difference remains constant at about 
20%, which is too large to use the approximation in cal- 
culations of the magnetization. In Session IIVI we show 
that this does not disqualify the chord susceptibility in 
calculation of the energy loss. 

In sections lnll and lTvl we present analytical solutions to 
the Debye relaxation equation enriched with the Larmor 
torque, 



dM M-M eq 

— - = hqjM x H q -, 

dt t 



(16) 



with M eq as given in eq. fTTj) . We will use these results 
to assess the effects of the torque on the steady-state so- 
lutions under linear and circular polarization of the ac 
magnetic field. As we take into account the curvature of 
the Langevin function, we can also discuss the implica- 
tions of the chord susceptibility. 



III. CIRCULAR POLARIZATION 

In this section we give the analytical solution of eq. (JTHJ) 
for a rotating magnetic field. The first term, represent- 
ing the Larmor torque, spoils the separation of Carte- 
sian components, which has enabled the derivation of 



the solution, eq. ([T2")) . found in Section [XT] for the De- 
bye relaxation equation. For a rotating magnetic field, 
H x = Hq cos(wi); H y = Hosm(ujt), the coupled equa- 
tions to be solved are as follows: 



dM x 

dt 
dM y 

dt 
dM z 

dt 



= -[j, jM z H sin(wt) - - (M x - M eq (H ) cos(wi)) ; 

T 

= hq^MzHq cos(wt) - - (M y - M eq (H ) sin(wt)) ; 

M, 

= Hol(M x H sin(wt) - M y H cos(wi)) -. (17) 



To handle the entanglement of the three components 
of M, it will prove to be convenient to apply the series 
of transformations that enabled us in a previous work 18 
to solve this set of equations for free precession, i.e. in 
the absence of relaxation (1/t = 0). Three transforma- 
tions create a coordinate system, which rotates as dic- 
tated by free Larmor precession in the rotating field. 
First, a rotation around the z axis by an angle lot makes 
the xy plane follow the magnetic field, then the z axis 
is turned by an angle into the direction of the to- 
tal angular velocity fl = lo + ljl, and finally a rotation 
around this new, z' axis by an angle fir drives the x'y' 
plane to rotate together with the magnetization vector. 
The description is reminiscent of an Euler transforma- 
tion and indeed the product of the three matrices repre- 
senting the rotations listed above, is of the form of the 
canonical Euler transformation^ with the replacements 
a f-> fit, /3 f-> — and 7 f-> — cut. As the axis of the 
Larmor precession is the magnetic field, which rotates in 
the plane perpendicular to the z axis, u) is perpendicular 
to wx, and SI = y 'ui 2 + u> 2 L . Also, it follows from this 
configuration that sin0 = uj^/VL and cosO = lu/VL. In 
what follows, the transformation matrix will be used in 
the form 



l 



' io cos tot cos Qt 



O — — lo cos Lot sin fit 

' ' lol cos tot 



f2 sin Lot sin f2i lo sin Lot cos fit — f2 cos Lot sin Sl< — lol cos Qt^ 
fl sin Lot cos Sit lo sin Lot sin fit + $1 cos Lot cos fit — lol sin fit 

lol sin Lot lo 



(18) 



To find the derivative of the transformed magnetiza- 
tion vector M', we need the derivative of the matrix O: 



equations for the transformed magnetization: 



dM' dOM dO^ dM 
= -=M + O- 



dt 



dt 



dt = dt 



(19) 



Substituting here from eq. lfT!))) . we find that the 



contribution of the Larmor torque to O 



^ cancels 



ilO 



ffM, as it should, leaving the following differential 



dM' 

dt 



OM + O 



'M eq {H )cos{Lotf 
M eq (H ) sin (Lot) 



(20) 



The first term on the right-hand side gives trivially 
M', the second one can be found applying eq. fTS)) . to find 
the differential equations for the three components of the 
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transformed magnetization vector. Ultimately, the ana- is 
lytical solution for the transformed magnetization vector 



M' x (t) = 
M' y {t) = 

Kit) = 



M' x (0) - M eq (H a ) 
M' V {Q) + M eq (H ) 



cos S 



y/l + (fir) 2 
to sin 5 

M' z (0)-M eq (Ho)^] exp(-</r). 



exp(-i/r)+M e9 (JJo) 
exp(-t/r) + M eq {H Q ) 



lo cos(fii — S) 

' ;! " ^7ttw ; 

lo sin(fii — 5) 



(21) 



Here 5 is defined by sine) = fir/ \/l + (fii~) 2 , cos<5 = 
+ (fir) 2 . Since the exponentially decaying terms 
are of no interest for applications on time scales exceeding 
r by several orders of magnitude ((/t> 1) and we seek 
the steady-state solution in the laboratory frame,we shall 
drop the exponential terms. The inverse transformation 
is easily carried out with the transposed of matrix (p~8|) . 
A compact form of the final result is 



E ~ A(tw) 



M(t) = M eq (H ) 



fi 2 1 + (fir) 



ujcos(ujt) + fi 2 r sin(wi) N 
to sin(wi) — fi 2 r cos(wi) 

(22) 
1 indi- 



Note that \M(t)\ = M eq (H )^l - ^ 

eating that the interplay of Larmor precession and relax- 
ation pushes the magnitude of the magnetization, which 
is time-independent, below the thermal equilibrium value 
in a static field Hq. The energy loss per cycle is easily 
calculated, 




FIG. 5: (Color online) The dependence of energy E ~ A = 
A(lot) on lot for different iol/io ratios as n — 1 + lol/lo. The 
curves reach maximum of 0.5n _1 at 1/n. Largest energy dis- 
sipation is achieved at ujl <C to. The maximum is shifted 
towards lower frequencies as n is increased, provided r is con- 
stant. 



f 0+27r/w ^jj 
E = -,/„ / M ■ — dt 



2irn M eq (H )Hc 



fir 



fi 1 + (fir) 2 ' 



(23) 



This Debye-type dependence on fir is familiar in low- 
field magnetic resonance^ and was shown by Garstens=I 
to be exact in the limit of vanishing static field. 

The basic functional dependence of E on Hq , to and r 
is dictated by the term j^jy ^ ■ Noting, that lo\ ~ _ff 2 , 
it has maximum either as a function of i?o, w, t or 
products lot and w^t. To facilitate comparison with 
Garstens — result for frequency dependence of E, we se- 
lect the dimensionless parameters lot and lol/lo and de- 
fine n = ^Jlo 2 l /lo 2 + 1 and A = |f r+ffep ■ Figure [5] shows 
A (which is proportional to £7) as a function of lot at var- 
ious value of n. The energy loss is seen to increase with 
lot up to a maximum at 1/n, the value of A at the max- 
ima being l/(2n). It is remarkable how the interplay of 



relaxation and Larmor precession suppresses the energy 
loss, in the case of circular polarization. 

The specific absorption rate (SAR) defined as energy 
loss per second and per kg, is the measure of energy losses 
relevant to applications: 



SAR = 



Elo _ fi Q M eq (H )H Lo 2 fir 
2^p~ ~ p ~fi 1 + (fir) 2 



(24) 



The function describing the dependence on lot and n 
is now 



B = 



1 L0 2 T 2 



Tl+n 2 (LOT) 2 ' 



(25) 



which grow monotonically with lot, saturating at B = 
7^2- Figure HJ shows that the saturation takes place at 
about lot = 1, meaning that increasing the frequency will 
not basically further improve the SAR. 
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SAR ~ B(™) 



IV. LINEAR POLARIZATION 




To find the equation of motion for the magneti- 
zation subjected to a magnetic field oscillating along 
the z axis, the following vectors have to be in- 
serted into eq.CEU): H(t) = (0, 0, H cos(ut)), M x 
H = j^(M y cos(uit),-M x cos(ujt),0) and M eq (H) = 

(0,0,4>M d £ cos(a;t))), where uj l = MoM^o- 

The matrix of the transformation that will eliminate the 
Larmor term of eq. (|16p depends on the sign of the gyro- 
magnetic ratio. For both cases we can write the matrix 

as 



FIG. 6: (Color online) The dependence of specific absorption 
rate SAR ~ B = B(ujt) on ujt for different ujl/uj ratios as 
n = 1 + ujl/uj- B saturates at ujl 3> 1, however most of 
change takes place up to ujt — 1. Largest specific absorption 
rate is achieved at <u (n — 1). 



a ± 



± sin g cos g ^ 
cos g ± sin g 
1, 



(26) 



where the upper sign equals that of 7 and g{t) — 
(ljl/lu) sin(wi). The equation of motion for M'(t) = 
Q ± M(t) is 



dM'(t) d£ d 



dt 



dt 



M + O 



dM dO 



± dt 



(1b M ± Mo|7l£ ± (M x H) 



: Q,(M - M eq {H)) 



(27) 



Taking into account that dg/dt = ujl cos(a;i) (irrespec- 
tive of the sign of ojl), the first two terms in eq. (|2~T|) 
cancel each other. The equation of motion for M'(t) in 
terms of the new coordinate system is easy to find, as 
Q_ ± M = M 1 , by definition, and M eq (H) has only a z 
component, which is evidently not affected by ± . What 
remains in the transformed frame is the set of equations 
of motion 



dM' x 

dt 
dM' v 

dt 
dM' z 
dt 



-M' ■ 



M' ■ 



ML - M S C 



' hqM^VHq cos(wi) 
V kf 



(28) 



The solution for the first two equations must be the 
familiar exponential relaxation to zero, 



K, v (t) = K, y (o)expHM- 



(29) 



Since ± separates the xy plane and the z 



axis, the 

transformation back to the laboratory frame can be done 
in two dimensions. As g vanishes at t = 0, we have 
M x (0) = -M'(Q) and M y (0) = M x (0), whence 



M(t) 



( ±M y (0) sin g + M x (0) cos g 
\ My (0) cos g qF M x (0) sin g 

cos[g T v(0)] 



exp(-t/r) 



M ± (0) 



T sin[ ff =f y(0)] 



exp(-t/r), (30) 



where 99(0) is the azimuth angle of Af(0) and Mj_(0) = 
We can determine the time depen- 



M|(0) + M2(0). 
dence of the angular velocity: 



<P 



dip 
~dt 



,-1 M y 



tan — ^ = tan I =p tan I — sin(wt) =p tp(0) 



M x 

=F— sin(wt)=F</?(0); 



=Fa;i cos(wt), 



(31) 



with the upper sign for 7 > and the lower sign for 7 < 0. 
Equation (|3Tj) describes a Larmor precession whose fre- 
quency is following the time dependence of the magnetic 
field (note that uj, is defined as the Larmor frequency 
at H — Ho, so that oj£cos(u;i) gives the Larmor fre- 
quency at H = H(t)). Whether the time dependence 
of the precession frequency is observable will depend on 
the magnitude of the projection of M on the xy plane, 
which, according to eq. (|30[) . decays exponentially, 



yjM*(t) + M*(t) = M ± (0) exp(-i/r). (32) 

If the relaxation time is very short (like in hyperther- 
mia, where lot w 10~ 4 , the magnetization will be fully 
aligned along the z axis before the magnetic field under- 
goes a substantial change, let alone a change of sign. Of 
course, to observe the time dependence given in eg. (1301) . 
we also need a Larmor frequency larger than the fre- 
quency of H, that is, lol > u>. This conditions is met in 
hyperthermia, at Hq = 200 kA/m the Larmor frequency 
is about 40 GHz, whereas tv is of the order of 100 kHz. 

The z component, being separated from those perpen- 
dicular to it, need not be put in a rotating frame. The 
search for an analytical solution to the last equation un- 
der (f2"5|) leads to a single (though not simple) integral if 
we substitute into eq. (j2"8")) the dM z /dt function extracted 
from 



d(^M,) = l et , et/r dM, 

dt t dt 



(33) 



HoHoMdV/kT <C 1 one can use the susceptibility 
instead of the Langevin function, i.e., the first term in 
the Taylor series of C instead of the function C. In fact, 
the integrals can be done for higher order terms also 
fRef. I22I, p. 228). As the Taylor series of cothx necessary 
here (Ref. I22I p. 42) is only valid for |x| < tt, the resulting 
expression for the time dependence of the equilibrium 
magnetization in terms of odd powers of cos(wi), 



£(( cos(w*)) = £ T^^WC cobM)) 2 "- 1 , (34) 

m=l ^ m '' 

where the -Bj™ are Bernoulli numbers as defined and 
listed in Ref. H^, p. 1040 and p. 1045, respectively. Substi- 
tuting the series (|34|) into eq. (|33l) . 



2m-l 



- (2m; 



gt/T COf ?™-^tdt. 

(35) 



To avoid the difficult integral, for £ = 



Now is straightforward to find a solution to the equa- 
tion of motion for M z which is analytical, albeit in the 
form of an infinite series : 



J 



Mz = 



2M 



s \ - 



— -E>2mC 



2 m- 



_1 E 

k=0 



(1/t) cos((2fc + l)ut) + (2k + l)wsin((2fc + l)ut) 
(m - 1 - fc)!(m + k)\ [(1/r) 2 + (2k + l)w 2 ] 



(36) 



There are three features of this solution that can be 
revealed without evaluating the series: (i) the exponen- 
tial factors have disappeared, meaning that the function 
describes a steady state, (ii) keeping only the m = 1 
term in the series the well-known result for the frequency- 
dependent Curie susceptibility [Ref. 11, eqs. @ and 
CEO]) ], M z = xoHo[cos(ujt) + wrsin(u>i)]/[l + (wt) 2 ] with 



Xo 



3kT 



(37) 



emerges and (iii) the double summation is essentially a 
Fourier series, of which only the k = sine terms are 
needed for the calculation of losses, because M z , mul- 
tiplied by —dH(t)/dt = Haujsm(ujt), will be integrated 
over a cycle. 

Keeping only the k — terms in the second sum in 
eq. (|36|) amounts to eliminating the higher harmonics in 
the time dependence of the magnetization, which are gen- 
erated due to the substantial nonlinearity of the Langevin 
function beyond x w 0.5. It follows then that for any 
value of Hq a susceptibility can be found, which will pro- 
vide the exact energy loss on the assumption that the 
magetization is 



M 



eff _ 



2M< 



E ( m l)2 
m=l v ' 



B2m( 



2 m -i cos(ujt) + lot sin(wi) 



1 + (wr) 2 



(38) 

Looking back at FigJ3](d), it becomes now clear that 
abrupt changes of the magnetization triggered by the flip- 
ping of the equilibrium magnetization will not influence 
the calculated energy loss, because they are represented 
by high-frequency terms in M z . 

The effective magnetization of eq. (38) may be looked 
upon as a fictitious magnetization, which, if realized, 
would give the correct losses without reproducing the 
correct fictitious hysteresis loop. Instead, the shape of 
the hysteresis loop and its frequency dependence are ex- 
actly the same as those one gets at low fields, where 
the static magnetization is proportional to the magnetic 
field. Equation (|3"5)) confirms then the insight that under- 
lies Rosensweig's introduction of the chord susceptibility 
[Refj 1 ^ eq.([2])] and enables the calculation of the "field 
dependence" of a fictitious susceptibility. In fact, the sus- 
ceptibility relevant to the energy loss, implicitly defined 
by Ml ff = xio SS H [cos(iot)+LOTsm(Lot)}/[l + (uJT) 2 } with 
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FIG. 7: (Color online) Convergence of the Taylor series ex- 
pansion of representing the summing term in eq. (|39p . 
Thin constant line is related to m ma x = 1 {xo), while thin 
curves are superposed functions up to mmax = 17. Red 
dashed line shows the exact xc = £(C)/C- Arrows point to- 
wards the direction of increasing m. 
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is not field dependent (that would contradict the defi- 
nition of the susceptibility as a limes at H — > 0), but 
amplitude dependent, through which is proportional 
to Hq. Note that B 2 — 1/6, so that xioss — Xo f° r C "C 1. 



In practice, the summation in eq. (|39[) will be limited to 
a finite index m max . The rapid increase of the Bernoulli 
numbers with m is a concern, but Figj7]shows that taking 
m max = 17, the convergence is reliable up tp C = 2.8. 
The converged values of £(£)/£ represent the amplitude 
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FIG. 9: (Color online) Specific absorption rate as a function 
of frequency to, parameterized by field strength Ho- Thick 
lines correspond to linear polarized cases, while dashed ones 
to circular polarized cases. Auxiliary thin line is drawn to ex- 
plore behavior around resonance lot = 1. Maximum reached 
for each circular polarized case when condition 10 > lol is 
fulfilled (Thin, dashed auxiliary lines). 



FIG. 8: (Color online) Energy loss as a function of frequency 
10, parameterized by field strength Ho- Thick lines corre- 
spond to linear polarized cases, while dashed ones to circular 
polarized cases. Auxiliary thin line is drawn to explore be- 
havior around resonance lot = 1. There is local maximum for 
each circular polarized case when condition lo = lol is fulfilled 
(Thin, dashed auxiliary lines). 



dependence of of the fictitious susceptibility Xloss- Note 
that while the deviation of Xchord from xo is increasing 
substantially, the difference between the exact xioss ancl 
Xchord is small and remains constant at about 5% even 
for 2 < ( < 2.8. Hence for ( > 1 calculating the energy 
loss using xo is not reliable, but close enough results can 
be obtained using the Xchord susceptibility. 



XI c 



t HoH M^ V4> 
kT 



E 



— ' (to!) 2 



B 2m C 2(m ~ 1 \ (39) 



Calculating the energy loss per cycle, 



E = - 



/'o 



tti+2-K /uj 



to 



„ r d H u 

M ■ ——dt 

dt 



r t a +2Tv/ul ITT 

no I Mff ■ ^dt 

'to 



dt 



KVoHoXiossiO 



TUl 



1 + (™) 2 



(40) 



r 



Comparing the lot dependence of E in both polar- 



ized cases, prompt can be realized that frequency de- 



10 



pendence of energy dissipation in circularly and linearly 
polarized cases can be equivalent by a factor of 2, if 
n = ^/ujl/uj + 1 = 1 and the linear approximation of 
M eq is taken for both polarized cases. 

At high driving field amplitudes (£ 3> 1) the energy 
absorption per cycle can be estimated as follows. For 
u) -C 1/r eq. (|28|) implies that M z lags behind the tar- 
get MsC(C cos(wi)) value by time r, which yields for the 
energy absorption E w A^lqMsHqlot. For ui ^ 1/r ac- 
cording to eq. ([28f the driving magnetic field causes only 
saw-tooth curve like oscillation of magnetization around 
the zero value with amplitude « Msit /2ujt. The energy 
absorption per cycle in this case is E rj AhqMsHq/lot. 

V. SUMMARY AND CONCLUSIONS 

Analytical solutions of Bloch-Bloembergen equation 
for both polarized cases are presented. They are in ex- 
cellent agreement with the numerical solutions. 

Differences in time-dependence of magnetization using 
Rosensweig's "chord" - susceptibility and Langevin func- 
tion are emphasized in Figs 13] and |4j where Hq is param- 
eterised. Formulae for energy dissipation and specific ab- 
sorption rate (SAR) are also provided. Short analysis is 
given for strong and weak field limit. It can be concluded 
that in hyperthermia region Rosensweig's suggestion for 
"chord" - susceptibility is reasonable. Linear approxima- 
tion of Langevin function, C{x) « x/3, can also be used 
provided the value of argument x is small enough. While 
x = 0.4 means 1% difference in linear approximation, at 
x = 1.2 already deviates with 10%. 

The frequency dependence of energy dissipation and 
SAR is also clarified, see Figs. [3] and [5] Possible maxi- 
mum of energy dissipation decreases with increasing ratio 
lul/lu and shifted to lower uj when relaxation time r is 
taken constant. SAR in the high lot limit shows satu- 
ration and its value is decreasing with increasing ratio 
ljl I lo . Equivalence of linear and circular polarized field 
can be stated with a factor of two, if the Larmor fre- 
quency is negligible small to external field frequency and 
linear approximation of Langevin function is taken for 
expressing equilibrium magnetization. In Fig(S] the re- 
sulting energy loss per cycle is depicted as a function of 
frequency for both polarized cases as a parameter of driv- 
ing field amplitude Hq where tq = 10 -9 s and Ms — 446 



kA/m with r — 5 nm nanoparticles. For small driving 
field amplitudes a common energy loss local maximum is 
observed. By increasing the driving field amplitude, the 
tent shape form of the energy loss curve is not changing. 
For the linearly polarized field it is proportionally shifted 
to higher energy loss values, while for the circularly po- 
larized case the energy absorption maximum is shifted 
towards higher frequencies. Notably, in the latter case 
at the given set of parameters the energy loss at lower 
frequencies just started to decrease with increasing field 
strength (see eq . (f2"B"l) ) . 

Figure [5] shows the same data in practical units of 
kWg -1 as defined for SAR. The presented results con- 
cerning the shift of the energy dissipation maximum 
and high frequency behaviour are in excellent agreement 
with the results derived from the Landau-Lifshitz-Gilbert 
equation^. 

Contrary the results obtained for low frequency region 
show different behaviour. This difference is outset by the 
possibility of changing the magnitude of magnetization. 
For example this produces resonance like E(Hq) depen- 
dence in Bloch-Bloembergen case, while the lack of this 
possibility results in Ho independent E(Hq) behaviour 
in the Landau-Lifshitz-Gilbert frame for circularly polar- 
ized field. The difference in results for linearly polarized 
field also originates from this possibility. 

At low frequency and relatively high driving field am- 
plitudes, linearly polarized field would yield more dissi- 
pation compared to the circular case when the relaxation 
processes characterized by the Bloch-Bloembergen equa- 
tion. Our final conclusion is alike as drawn from analyz- 
ing Landau-Lifshitz-Gilbert equation: when the relax- 
ation process is characterized by the Bloch-Bloembergen 
equation, the technical complications of generating a cir- 
cularly polarized field in difficult geometry need not be 
considered. 
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